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ABSTRACT 

We study the dust depletion pattern in eight well separated components of the 
Zabs = 1-973, log N{R i) = 20.83, damped Lyman-a system toward Q 0013-004, 
four of which have detectable H2 absorption. The apparent correlation between the 
abundance ratios [Fe/S] and [Si/S] in the components indicates that the abundance 
pattern is indeed due to dust-depletion. In particular, we find evidence for deple- 
tion similar to what is observed in cold gas of the Galactic disk ([Fe/Zn] = —1.59, 
[Fe/S] = -1.74, [Zn/S] = -0.15, [Si/S] = -0.85) in one of the weakest components 
in which molecular hydrogen is detected with log 7V(H2) ~ 16.5. This is the first time 
such a large depletion is seen in a DLA system. Extinction due to this component is 
negligible owing to small total Hi column density, log iV(Hi ) ^ 19.4. This observation 
supports the possibility that current samples of DLA systems might be biased against 
the presence of cold and dusty gas along the line of sight. 

The overall metallicities of this peculiar DLA system in which O i and C ii are spread 
over ~1050 km s^^ are [P/H] = -0.64, [Zn/H] = -0.74 and [S/H] = -0.82 relative to 
solar. The sub-DLA system at Zabs = 1.96753 has [P/H] > 0.06, [Zn/H] > -0.02 and 
[S/H] > —0.18. The overall molecular fraction is in the range —2.7 < log / < —0.6, 
which is the highest value found for DLA systems. H2 is detected in four com- 
ponents at —615, —480, and 85 km s^^ relative to the strongest component at 
Zabs = 1.97296. CO is not detected (log A^(CO)/iV(H i) < -8) and HD could be 
present at Zabs — 1.97380. 

We show that the presence of H2 is closely related to the physical conditions of the gas: 
high particle density together with low temperature. The observed excitation of high J 
H2 levels and the molecular fraction show large variations from one component to the 
other suggesting that the UV radiation field is highly inhomogeneous throughout the 
system. Gas pressure, estimated from C i absorptions, is larger than what is observed 
in the ISM of our Galaxy. This, together with the complex kinematics, suggests that 
part of the gas is subject to high compression due to either collapse, merging and/or 
supernovae explosions. This is probably a consequence of intense star-formation ac- 
tivity in the vicinity of the absorbing gas. 

Key words: Cosmology: observations - Galaxies: halos - Galaxies: ISM - Quasars: 
absorption lines - Quasars: individual: Q 0013—004 



1 INTRODUCTION 



* Based on observations carried out at the European Southern 
Observatory (ESO) under prog. ID No. 66.A-0624 and 267.A- 
5714 with UVES the echelle spectrograph mounted on the Very 
Large Telescope (VLT) at Cerro Paranal Observatory in Chile. 



The amount of dust present at high redshift has impor- 
tant consequences on the physics of the gas. In addition, 
dust directly affects our view of the high redshift uni- 
verse through extinction. Therefore, the presence of dust 
in damped Lyman-a (hereafter DLA) systems, that contain 
most of the neutral hydrogen in the universe, can have sig- 
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nificant consequences. Although the presence of dust in DLA 
systems has been claimed very early (Pei et al. 1991), the is- 
sue has remained controversial. Indeed, Lu et al. (1996) have 
questioned the idea that the overabundance of Zn compared 
to Cr or Fe observed in DLA systems (e.g. Pettini et al. 1997) 
is due to selective depletion onto dust-grains and have ar- 
gued that the overall abundance pattern observed in DLA 
systems is indicative of Type II supernovae enrichment in- 
stead. In recent years several studies have shown that both 
effects, dust-depletion and peculiar nucleosynthesis history, 
should be invoked to explain the abundance pattern (Vladilo 
1998, Prochaska & Wolfe 1999, Ledoux et al. 2001a). How- 
ever, the lack of statistics and the wide variety of objects 
that can give rise to DLA systems, namely dwarf galax- 
ies (Centurion et al. 2000), large disks (Prochaska & Wolfe 
1997, Hon et al. 2001), galactic building blobs (Haehnelt et 
al. 1998, Ledoux et al. 1998) etc., with, for each of these ob- 
jects, its own history, prevent us to have a clear picture of the 
nature of DLA systems. Nevertheless, all studies conclude 
that the dust content of DLA systems is small. However, it 
is possible that the current sample of DLA systems is biased 
against high-metallicity and dusty systems. Indeed, Boisse 
et al. (1998) have noticed that there is a lack of systems 
with large A'^(H l) and large metallicity. Very recent investi- 
gation of an homogeneous sample of radio-selected quasars 
shows that the dust-induced bias cannot lead to underesti- 
mate the H I mass in DLA systems by a large factor (Ellison 
et al. 2001). However even a factor of two could change our 
understanding of DLA systems. 

An obvious way to search for DLA systems with large 
amount of dust is to select those where molecules are de- 
tected as these molecules form predominantly at the surface 
of dust grains. However, it may not be so simple as it has 
been shown that the presence of H2 is not only related to the 
dust-to-metal ratio but is mostly dependent on the physical 
conditions of the gas. First of all, H2 is detected when the 
particle density is large (Petitjean et al. 2000, Ledoux et 
al. 2001b). In any case, the system at Zabs = 1.973 toward 
Q 0013—004 is a good target as very strong molecular ab- 
sorption lines have been identified by Ge & Bechtold (1997, 
see also Ge et al. 1997). 



2 OBSERVATIONS 

The Ultra-violet and Visible Echelle 

Spectrograph (D'Odorico et al. 2000) mounted on the ESO 
Kueyen 8.2 m telescope at the Paranal observatory has been 
used on October 21—23, 2000, to obtain high-spectral reso- 
lution spectra of Q 0013—004. Additional observations were 
performed in service mode on September 10—11 and 13—18 
2001. The slit width was 1 arcsec (the seeing FWHM was 
most of the time better than 0.8 arcsec) and the CCDs were 
binned 2x2. The resulting resolution was ~45000. The to- 
tal exposure time ~22 hours was split into 1 h exposures. 
Two dichroic settings were used to cover the complete wave- 
length range. The data were reduced using MIDAS, the ESO 
data reduction package, and the UVES pipeline in an inter- 
active mode. The main characteristics of the pipeline is to 
perform a precise inter-order background subtraction, espe- 
cially for master flat-fields, and an optimal extraction of the 
object signal rejecting cosmic ray impacts and subtracting 



Table 1. Fit results for H2 lines 
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the sky at the same time. The data reduction is checked step 
by step. Wavelengths were corrected to vacuum-heliocentric 
values and individual ID spectra were combined together. 
The resulting S/N ratio per pixel in the final spectrum is of 
the order of 35 at ~3500 A and 80 at ~5000 A. 



3 THE DLA SYSTEM AT Zabs = 1.973 
3.1 Molecular hydrogen 

The presence of H2 in the z^hs = 1.973 system to- 
ward Q 0013—004 was first reported by Ge & Bechtold 
(1997). Based on an intermediate resolution spectrum and 
a single component curve of growth analysis, they de- 
rived a very high value of the molecular fraction, f — 2 
N{U2)/[N{U i)+2iV(H2)] = 0.22 ±0.05. In the UVES data, 
H2 is detected in four distinct main components spread 
over ~700 km s"\ at Zabs= 1.96685 (~ -615 km s"^), 
1.96822 (~ -480 km s'^), 1.97296 (0 km s'^) and 1.97380 
(~ -1-85 km s^i), marked with vertical dotted lines in Fig. |l| 
For simplicity we name these components a, b, c and d 
respectively. Four additional strong metal components at 

Zabs= 1.96729 ( 572 km s^^), 1.97023 ( 275 km s'^), 

1.97098 (~ -200 km s'^) and 1.97138 (~ -159 km s'^), 
with highly saturated O I, Si 11 and Mg 11 lines, but clearly 
identlGed by the weaker S 11 lines, are marked with vertical 
dashed lines in Fig. ^(see also Fig. ^and Fig. 4). These addi- 
tional components, in which we measure a conservative up- 
per limit, A'^(H2)^ 10^* cm^ , probably have a non-negligible 
contribution to the total H I column density. 
In the case of components c and d the absorption lines due 
to J 3 rotational levels are strongly saturated and possi- 
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Figure 1. H2 absorption profiles for transitions from difl^erent levels, v = — (left panels) and 1 — (right panels), given on a velocity 
scale with origin at 2abs= 1.97296. The vertical dotted lines indicate the positions of the four H2 components discussed in the text. The 
location of components with strong O I, Si II, Mg II and S II absorption lines but without detectable absorption due to H2 is indicated 
by dashed vertical lines (see also Fig. 2). 



bly blended with other lines. This implies that we can only 
derive a lower limit on the H2 column density from the con- 
dition that the lines must be saturated and an upper limit 
from the absence of damping wings. Results are given in 
Table |l| together with the number of unblended lines, nuncs, 
that are used to derive the column densities. 
Components a and b are narrow and mostly unblended. As 
the lines are moderately saturated it is possible to derive 
H2 column densities for various J levels (see Table |l|). The 
total H2 column density in these components is 2.39±0.22 x 
10^*^ and 3.46±0.40 x 10^® cm"^ respectively. The excitation 
temperatures Tqi = SOOl'lgg^ and K could reflect the 

kinetic temperature of the gas (e.g. Abgrall et al. 1992). 
In the case of component a the excitation temperatures for 
the ortho and para states (T02 = 776II44 K and Tn = 
40lt^5 K) are diflterent (although errors are large) . This hap- 
pens also in component b. Note however that the estimated 
column densities of J = and ,1 = 2 levels are based on 
a single unblended line. Moreover, J = transitions from 
this latter component are blended with J = 1 transitions 
at Zaba=l. 96685 (see Fig. |). The fact that T02 and Tn are 
different could indicate that excitation is not dominated by 
collisions and other processes such as UV pumping are im- 
portant (see Srianand & Petitjean 1998). As emphasized by 
Srianand et al. (2000), populating the J = 4 level after for- 
mation of a molecule on a dust grain is probably negligible 



in this gas. The photo-absorption rate in the Lyman and 
Werner bands is therefore /3o ^ 10~^ and 9x10"^^° s"'^, re- 
spectively, for a and b, based on the J = 4 column densi- 
ties. This is consistent with what is measured in the Galactic 
ISM. 

Absorption due to H2 is detected up to J = 5 in component 
c. The lines from the J = 4 and 5 levels are narrow and 
unsaturated and we used a single-component fit to derive 
the column densities. Neglecting the formation pumping to 
J = 4 and using the lower limit on the J = column density, 
we obtain (3o ^ 3.0 x lO'^^s"^. This is at least a factor 3 
smaller than that derived for component b. 
It is clear from Fig. |l| that the J ^ 4 lines are either absent or 
weak in component d. We determine the column density for 
J = 3 from the H2-LO-R3 line that is not saturated. Absence 
of absorption from high J levels, despite strong saturation 
of the transitions from J ^ 2 levels, suggests that the rate 
of UV pumping to high J levels is low in this component. 
We derive /3o ^ 1.1 x 10~^''s~^ which is at least a factor 8 
smaller than what is derived for component b. 
All this indicates that there are large fluctuations of the 
ambiant UV flux inside the system. 
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Figure 2. Metal line profiles on a velocity scale with origin at rcdshift z^bs = 1.97296. The vertical dotted and short-dashed lines are 
as in Fig. ^. The vertical long-dashed lines indicate narrow components seen only in low-ionization species. For doublets, the profile of 
the weakest line is overplotted as a dashed line. 



3.2 Atomic hydrogen 

The overall structure of the strong metal lines is illustrated 
in Fig. ^ The location of the strongest components with and 
without detected H2 molecules is indicated by, respectively, 
dotted and short-dashed lines as in Fig. 0. The long-dashed 
lines show components with absorption due to low-ionization 
species only. It can be seen that the low-ionization absorp- 
tions (from O I, C II, Si 11, Mg 11) span about 1000 km s"^ and 
are structured in two main strong systems at Zabs ~ 1.9676 
and 1.9733 separated by 575 km s~^. Note that most of the 
Si IV and C iv absorption is located inbetween these two 
strong low-ionization systems. 

A Voigt profile centered at the mean position of the two 
strongest H2 components (c and d), Zahs = 1.9733, gives a 
good fit to the red wing of the Lyman-a absorption with 
a total column density, iV(H l) = 6.7x10^" cm"^. This is 
consistent with the value derived by Ge & Bechtold (1997). 
There is however an excess of absorption in the blue wing 
that we attribute to H2 components a and b that are blue- 
shifted by 615 and 480 km s^^ relative to the strongest H2 
component. We can therefore estimate A''(H l) associated 



with these two blueshifted components. In Fig. |3|, we show 
the fits to the damped absorption with and without the ad- 
ditional absorption (log A''(H l) = 19.4) that can be accomo- 
dated at 2abs ~ 1.96753. This implies that less than 10% of 
the Hi gas is associated with H2 components a and b. This 
gives a lower limit on the molecular fraction / > 2 x 10~^ 
and 4x 10""^ for components a and b respectively. These lim- 
its are higher than what has been observed in other DLA 
systems (e.g. Petitjean et al. 2000). Such a large value for 
/ in a low Hi column density cloud is at odd with what is 
observed in the ISM of our Galaxy where log / < —4 for 
log iV(H l) < 20 (Savage et al. 1977). This implies high for- 
mation rate, and therefore large amount of dust, high den- 
sity, low temperature, and low destruction rate, therefore 
shielding from the UV field. 

We derive for the whole system 6x10^^ < log N{R2) < 
1 X 10'^'' cm~^ . Therefore, the mean fraction of molecular hy- 
drogen is in the range 0.002—0.25. This is much larger than 
what is seen in other DLA systems (see Petitjean et al. 2000; 
Ledoux et al. 2001b). 
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Figure 3. Fits to the damped Lya absorption profile at 
2abs = 1-9733. Fits with (solid line) and without (dotted line) 
the maximum column density (log N{R l) = 19.4) present in the 
H2 components a and b at a mean redshift z^bs = 1.96753, com- 
patible with the asymmetry of the profile. 



3.3 Other molecules 

If the molecular fraction is close to the upper limit 
we derived above then one can expect to detect other 
molecular species such as CO and HD (see Varshalovich 
et al. 2001). We do not detect molecular CO down to 
iOobs(A1088) < 7 mA. This means N{CO) < 2x10^^ cm"^ in 
components c and d, using the oscillator strength by Mor- 
ton (1975). Therefore log iV(CO)/iV(H i) < --8.5. Note that 
in our Galaxy, when log A'^(H2) > 19, CO is readily detected 
with log iV(CO) ^ 13 (Federman et al. 1980). 
It is difficult to trace HD in this system because blending of 
lines is severe for all bands. At Zabs = 1.97296 (component 
c), absorption features are present at the expected position 
of HD-LO-RO and HD-Ll-RO but they are stronger than a 
possible HD-L3-R0 absorption {woha ~ 30 mA). Note that 
consistent features are seen at 2:abs = 1.97278, especially at 
the expected position of HD-L3-R0 but HD-LO-RO is not 
visible down to tUobs < 40 mA. The large equivalent width 
limit is a consequence of the line falling in the wing of an- 
other stronger absorption. At Zabs = 1.97380 (corresponding 
to component d), consistent absorption features are present 
at the expected position of HD-LO-RO, HD-L3-R0 and HD- 
L4-R0. HD-Ll-RO is possibly present but blended. The de- 
tection should be considered as tentative as HD-LO-RO and 
HD-L3-R0 are affected by noise. The strongest constraint 
on this component is given by Wobs (HD-LO-RO) ~ 10 mA. 
Using A = 1105.845 A and / = 7.6x10"* we derive 
A''(HD) ~ 4x10^* cm~^. Using the limits on H2 indicated in 
Table 1, we derive 1.0x10"^ < iV(HD)/iV(H2) < 2xl0~^ 



4 METALLICITY AND DEPLETION 
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Figure 4. Profiles of various metal lines detected in the DLA 
system plotted on a velocity scale with origin at Zabs = 1.97296. 
Solid lines give the best Voigt profile fits. Vertical dotted lines 
indicate the positions of subcomponents. 



4.1 The system as a whole 

The absorption lines from heavy elements show a complex 
multiple component structure (see Fig. ^ and only a few of 
the strong velocity components have detectable associated 
H2 absorption. 

The column densities of different species, integrated over, 
respectively, the DLA system at Zabs = 1.9733 and the 
sub-DLA system at Zabs = 1.96753 are given in Table ^ 
For the sub-DLA system, we have considered that most 
of the neutral hydrogen is associated with the two molec- 
ular components a and b. The column densities indicated 
in Table ^ are therefore the sum of the column densities 
measured in these two components. Note that the measure- 
ments are obtained using transitions with little saturation 
(see Fig. 4). We use in the following the standard notation 
[X/H] = log[Z(X)]-log[ZQ(X)] with Z{X) the metallicity of 
species X. Solar metallicities are from Savage & Sembach 
(1996). 

Sulphur and zinc are not much depleted onto dust-grains 
in the ISM of our Galaxy and therefore are good indica- 
tors of metallicity. They have consistent metallicities, about 
0.15 solar in the DLA system and 0.7 solar in the sub-DLA 
system. Phosphorus is slightly enhanced (~ 1.5 times) com- 
pared to Sulphur and Zinc in both cases. Iron, Chromium, 
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Table 3. Column densities and metallicities relative to solar in 
the Zabs = 1.96822 component 



■100 100 

velocity (kms"') 



Figure 5. Profiles of various metal lines detected in the H2 com- 
ponent at Zabs= 1.96822 (component b). Solid lines give the best 
Voigt profile fits. 



Table 2. Total column densities and metallicities relative to solar 



Ion 


N 




[X/S]''''^ 






(cm-2) 


DLA 


Cold 


Warm 


Halo 




DLA at Zabs = 1.9733'= 






H I 


6.70±0.70 X 10^" 










Si 11 


2.70±0.20 X W-' 


-0.13 


-1.3 


-0.4 


-0.3 


S 11 


1.90±0.10 X 10^5 


-0.00 




+0.0 


+0.0 


P 11 


5.71±0.77 X 10^3 


+0.18 


-0.5 


-0.2 


-0.1 


Zn 11 


5.47±0.50 X 10^2 


+0.08 


-0.4 


-0.2 


-0.1 


Cr 11 


9.86±0.30 X 10^2 


-0.70 


-2.1 


-1.2 


-0.6 


Fe 11 


6.41±0.46 X 10^* 


-0.71 


-2.2 


-1.4 


-0.6 


Ni II 


1.44±0.60 X 10^3 


-1.10 


-2.2 


-1.4 


-0.6 


Mn II 


2.31±0.50 X 10^2 


-1.18 


-1.5 


-1.0 


-0.7 




Sub-DLA at 




1.96753'' 






H I 


<2.70±0.20 X I019 










Si 11 


2.65±0.31 X 10^'' 


-0.38 








S 11 


3.32±0.29 X 10^* 


-0.00 








P 11 


1.16±0.17 X IO13 


+0.24 








Zn II 


1.16±0.12 X 10^2 


+0.16 








Fe II 


5.43±0.35 X 10^3 


-1.03 









^'[X/H] =log[Z(X)]-log[Z0(X)] 
'^Galactic values are from Welty et al. (1999) 
c[S/H] = -0.81 in the DLA system 
'^Components a and b; [S/H] > -0.18 



Species 


N 


[X/S] 


[X/S] 




(cm-2) 


Component 


Cold ISM 


Mgl 


3.58+0.20 X 10" 






Mg II 


1.24+0.10 X 10^^ 


-1.33 


-1.2 


Fe 11 


4.07+0.23 X 10^2 


-1.74 


-2.2 


Zn 11 


2.23+0.10 X 10" 


-0.15 


-0.4 


Si 11 


3.51+0.38 X 10^=* 


-0.85 


-1.3 


S II 


1.30+0.10 X 10'"' 


+0.00 


0.0 


I 


> 6.0 X 1014 






N I 


^ 6.4 X IQi-i 






P II 


^ 3.0 X 10^2 


<0.06 




Al 11 


3.93+0.57 X 10" 


-1.73 


-2.4 



Nickel and Manganese are underabundant with respect to 
Zinc probably because they are depleted onto dust-grains. 
If we compare the measurements to what is seen in our 
Galaxy, it is interesting to note that the abundance pat- 
tern seen in the DLA system is close to that of warm halo 
gas whereas the pattern in the sub-DLA system is close to 
that of warm neutral gas in the disk (see however next Sec- 
tion). This is consistent with the conclusion by Petitjean et 
al. (2000) that the bulk of the DLA population is drawn 
from warm gas. 



4.2 Dust depletion in the z^bs = 1.96822 
component 

In contrast to what is seen for the system as a whole in which 
absorption by Fe 11 is strong compared to absorption by 
Zn II, it is stricking to note that in the Zabs = 1.96822 com- 
ponent (H2 component b), the absorption in the strongest 
iron transition Fe llA2382 has about the same strength as 
that in Zn llA2026 (see Fig. ||). In fact, all transitions from 
refractory elements are weak. Note that Si ill absorption is 
clearly absent (see Fig. indicating negligible ionization 
correction. The results of Voigt profile fits to the available 
metal lines are given in Table H Abundances of Zn and S are 
similar. As can be seen from Table |^, Fe and Al are depleted 
compared to Zn by about two orders of magnitude; Mg and 
Si are depleted by an order of magnitude compared to S. 
This is the first time that such a large depletion, consistent 
with that observed in the cold interstellar medium of the 
Galactic disk, is observed in a DLA system. 



4.3 Variations of dust-depletion within the system 

We investigate here the variations of the abundance ratios 
among different well separated velocity components where 
ionization correction is negligible. This is the case for the 
components where H2 is detected. However, heavy element 
absorption from component a is blended with absorption 
due to a nearby high- ionization feature (see Fig Thus 
we do not consider this component for the analysis. In addi- 
tion, we consider the two components at ~ —153 and —195 
km s"'^ that account for 20 to 30% of the observed heavy 
element column densities (see Fig. 4). 
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Figure 6. Depletion of Fe (stars) and Si (squares) with respect 
to S observed in different velocity components. The zero of the 
velocity scale is at z^bs = 1.97296. The values averaged over the 
components are indicated by horizontal dashed lines. 



We note that the [Zn/S] ratio is very close to solar (uncer- 
tainties are typically of the order of 0.1 dex) in all compo- 
nents. Si and Fe abundances with respect to S are plotted 
in Fig. ^ for different velocity components. Depletion of Si 
and Fe is apparent. All components, except the very spe- 
cial component b (at —480 km s~^ see previous Section), 
have similar abundance ratios (and therefore depletion fac- 
tors), irrespective of whether H2 molecules are present or 
not. This has already been noted in the ^abs = 1.962 DLA 
system toward Q 0551—366 (Ledoux et al. 2001b). 
In Fig. 1^ [Fe/S] is plotted as a function of [Si/S] as observed 
in the sub-components. It is apparent from the figure that 
there is a strong correlation between the two quantities. Note 
that the correlation is present (albeit with lower significance 
level) even if we do not consider the Zabs = 1.96822 compo- 
nent. This clearly demonstrates the existence of differential 
dust-depletion in a single DLA system. The depletion factors 
for three gaseous components of our Galaxy are overplotted 
on Fig. 1^: cold and warm ISM and halo gas (Welty et al. 
1999). It is clear that most of the velocity components have 
depletion patterns very close to that seen in warm and halo 
gas in our Galaxy. Only the Zabs = 1.96822 component has 
a pattern similar to that of the cold ISM phase. 
It is important to note that the lowest depletion factor is 
recorded in component d (at -|-85 km s~^) that has high H2 
content. In addition, components with similar dust-depletion 
factors can have very different H2 contents. This means that 
the presence of dust, though important to form H2 molecules 
at low temperature, is only one of the parameters, together 
with high particle density, low temperature and low UV- 
flux, that decide what is the molecular fraction in the gas. 



Figure 7. [Fe/S] vs [Si/S] for different velocity components. The 
corresponding values measured in different gas components in our 
Galaxy are indicated by stars. 



5 PHYSICAL CONDITIONS IN THE GAS 



C I is usually seen in components where H2 is detected (see 
however Srianand & Petitjean 1998) and it is possible to 
derive constraints on the particle density of the gas from 
the relative populations of the different levels of the C I 
ground term. 

Absorption lines produced by C I and C I* are clearly seen in 
the four components of the DLA system toward Q 0013—004 
in which H2 is detected. However, C I** is seen only in 
two components at ^abs = 1.96822 and 1.97296. In addi- 
tion, there are three components that show C I absorption 
without detectable H2 (^5 lO^'^cm"^, see also Srianand & 
Petitjean 2001). We perform simultaneous Voigt profile fits 
to the Ar = 1560 and 1656 A C l, C l* and C l** lines. The 
fit results are shown in Fig ^ and summarized in Table ^. 
Upper limits on the Cosmic Microwave Background Radi- 
ation (hereafter CMBR) temperature, Tcmbr, are derived 
assuming that CMBR pumping is the only excitation mech- 
anism. Results, given in column #7 of Table § are con- 
sistently larger than the value expected from standard hot 
big-bang model 8 K, see Srianand et al. 2000). This is 
because other excitation processes are at play and in particu- 
lar excitation by collisions. We estimate the particle density 
in the cloud assuming, Tcmbr = 8 K, Tkin = 100 K and 
neglecting fluorescence. The range given in column #8 of 
Table ^ are obtained using 2a uncertainties in the column 
density estimate. 

As the ionization potential of C I is similar to the energy 
of photons in the molecular bands, it is expected that C I 
and H2 follow each other. It is therefore interesting to note 
that A'^(C l) in component d is about a factor of 2 smaller 
than in component b whereas N{Y{2) is at least a factor of 
ten larger. This is probably a consequence of much larger 
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Table 4. Results of Voigt profile fitting to the C I absorption lines 



^abs 



N{C i) 
(10" cm-2) 



Af(C I*) 



(IQi 



N{C I**)" 



(10l2 cm-2) (km s-l) (10^ 



N{S II) 



TcMBR 

(K) 



(cm~^) 



96679 
96691 
96706 



11.6±1.00 
7.29±0.45 
2.06±0.34 



Component a 
5.23±0.58 
3.56±0.59 
1.80±0.45 



Af(H2)= 



3.10±0.57 
3.01±0.20 
4.93±1.44 



1.87±0.20 



^13.8 
^14.5 



20-70 
20-80 



96822 



10.9±0.03 



Component b 
12.6±0.06 



Ar(H2)=3.46±0.40 x lO^'^ cm" 



5.83±0.30 



4.30±0.20 



1.44±0.06 



^20.0 



170-200 



Component c - Af(H2)=0.03-6 X 10^^ cm'^ 

97280 11.9±0.03 5.89±0.35 .... 5.98±0.28 1.42±0.18 ^13.0 40-60 

97296 27.3±0.05 16.0±0.05 3.39±0.32 6.23±0.16 5.68±0.34 ^13.5 50-65 

97316 4.56±0.34 2.13±0.46 .... 7.00±0.77 0.82±0.10 ^15.0 10-85 



97365 
97382 
97399 
97417 



2.10±0.24 
2.52±0.23 
0.81±0.23 
0.31±0.16 



Component d 
^1.30 
1.10±0.31 



Ar(H2)=0.02-4 x lO^^ cm 
5.66±0.41 
4.52±0.42 
7.30±1.02 
0.37±0.27 



0.73±0.13 
1.32±0.15 



^15.0 
^15.5 



< 85 

< 95 



Components with N{U2)^ IQi^cm^^ 
96737 8.15±0.63 5.91±0.79 .... 14.41±1.08 

96763 7.08±0.83 4.76±0.99 .... 27.14±3.38 

97109 1.90±0.12 ^1.30 .... 10.05±0.91 

97144 9.10±0.62 8.05±0.60 .... 13.44±0.96 



1.72±0.19 
1.33±0.20 



< 17.5 
^ 19.3 
^ 21.0 
^18.5 



40-135 
20-180 
^ 200 
65-160 



1. 
1. 
1. 

1. 

1. 
1. 
1. 

1. 
1. 
1. 
1. 



1. 
1. 
1. 
1. 

°- The 2a upper limit is 1.5x10^^ cm"^. 

A'^(H2)/A''(H l) ratio in component d. This conclusion is fa- 
vored by the indication that the H I column density is prob- 
ably of the same order of magnitude in both components as 
indicated by the S ii column densities (see Table m. 
C I, C I* and C I** are clearly detected in component b (see 
Fig ^). The excitation temperature of the J = 1 H2 level is in 
the range 60-80 K. Assuming that the kinetic temperature 
has the same value and Tcmbr = 8 K, we derive that the 
particle density is in the range nu ~ 170-200 cm~'^. 
It is interesting to note that the components at 
Zabs = 1.96737, and 1.96763 have a large pressure. Consid- 
ering the kinetic temperature to be 100 K implies P/k — 
4000 - 13500 and 2000 - 18000 cm"^K respectively. How- 
ever, these components belong to the strong saturated ab- 
sorption component at ~— 600 km s~^ (see Fig. 2) with 
large O I and Si ill absorptions and no H2. It is most 
likely that the temperature of the gas is higher than what 
we use here (see Petitjean et al. 2000). For T = 1000 
K, we derive n — 10—75 cm~^ and 28—54 cm""*, and 
P/k = 18,000 - 61,000 and 10,000 - 75,000 cm'^K re- 
spectively for the two components. 

Such high pressures are seen in about one third of the C I 
gas in our galaxy (see Jenkins & Tripp 2001). This probably 
indicates that part of the gas in the system under consider- 
ation is subject to high compression due to either collapse, 
merging and/or supernovae explosions. The overall velocity 
structure of the system at z^hs = 1.9676 with the presence 
of molecular gas at -1-100 and —150 km s~^ from the center 
of the high-ionization gas strongly suggests the expansion of 
a spherical shell. In any case, intense star-formation activity 
is occuring in the vicinity of the gas. 



6 CONCLUSION 

The DLA system at Zabs = 1-973 toward Q 0013-004 is pe- 
culiar in several aspects. Absorption lines from metal species 
are spread over about 1000 km s~^; in particular O I and 
C II span 1050 km s~^. The velocity structure indicates the 
presence of two main sub-systems centered at z ~ 1.973 
and 1.9674, separated by ~ 550 km s~^ with, respectively, 
log N{U 1) = 20.83 and <19.4 and [S/H] = -0.81 and 
> -0.18. 

The low-ionization gas is conspicuous in the system. There 
is clear evidence that all species are at most twice ion- 
ized in the Zabs ~ 1.96822 component. This means that 
if photo-ionization dominates, there is probably very few 
photons with energy larger than 20 eV. More generally, 
N{X+)/N{X^+) > 1 over the entire z ~ 1.973 system. This 
ionization state could reveal gas ionized by slow shocks. This 
idea is reinforced by the high pressure measured in a few 
components from C I absorptions. 

H2 is detected in four main components, two very strong 
components (log A*'(H2) > 17) in the z ~ 1.973 sys- 
tem and two weaker components (log N{H2) ~ 16) in 
the z ~ 1.9674 system. The total column density is 
17.8 < log A'^(H2) < 20.0 and therefore the mean molec- 
ular fraction, / = 2iV(H2)/(2iV(H2) + iV(H l)), is in the 
range —2.7 < log / < —0.6 which is the highest molecular 
fraction observed in DLA systems. 

The analysis of the A'^(C l*)/A''(C l) ratio in different com- 
ponents indicates that whenever II2 is detected, the particle 
density is high {nu > 30 cm~^). High density is also found 
for components without any H2 absorption detected. There 
is a hint for the depletion of metals within the components 



Dust Depletion in Damped Lyman-a systems 9 




-600 -400 -200 200 



Figure 8. Velocity plot showing the absorption profiles in transi- 
tions from different levels of the C I ground-term. The continuous 
curve gives the best fitted profile for multiple components. The 
zero of the velocity scale is at z^bs = 1-97296. 

to be correlated to the A''(C l*)/A''(C l) ratio. This suggests 
that depletion onto dust-grains could be larger for denser 
gas. 

The component at Zabs = 1.96822 shows evidence for 
large depletion of iron and silicon relative to sulfur and 
zinc ([Fe/Zn] = -1.59, [Fe/S] = -1.74, [Zn/S] = -0.15, 
[Si/S] = —0.85) similar to what is observed in cold gas 
of the Galactic disk. The corresponding dust-extinction 
is small in this case because, although H2 is detected 
with log A''(H2) ~ 16.5, the H l column density is small, 
log A'^(H l) < 19.4, in the component. This is direct evi- 
dence for a considerable fraction of heavy elements being 
locked into dust-grains, and, as a consequence, this supports 
the idea that the current sample of DLA systems might be 
biased against the presence of cold and dusty gas along the 
line of sight. Note that the rest of the gas shows a depletion 
pattern close to that of warm gas in the Galactic halo. 
The overall kinematics of the system with the two strong 
sub-systems separated by 550 km s~^ suggests that the 
line of sight is passing through one or several objects in 
strong interaction. The velocity structure of the subsystem 
at Zabs = 1.9676 with the presence of molecular gas at +100 
and —150 km s~^ from the center of the high-ionization ab- 
sorptions strongly suggests the expansion of a spherical shell. 
All this, together with the strong inhomogeneity of the UV 



field, the high pressure in the gas and the high metallicities 
strongly suggests that intense star-formation activity is oc- 
curing in the vicinity of the system which should be revealed 
by deep imaging in the field. 
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